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CEPSTRAL DOMAIN PULSE OXIMETRY 

RELATED APPLICATION INFORMATION 
This application is a continuation of and claims priority from United States 
5 Application Serial No. 10/371,658 entitled "CEPSTRAL DOMAIN PULSE 

OXIMETRY" filed on February 21, 2003, which claims priority from United States 
Provisional Application Serial No. 60/359,018 entitled "CEPSTRAL DOMAIN PULSE 
OXIMETRY" filed on February 22, 2002, the entire disclosures of which are 
incorporated herein. 

10 

FIELD OF THE INVENTION 
The present invention relates generally to pulse oximetry, and more particularly to 
pulse rate and blood analyte level estimation using cepstral domain processing of 
plethysmographic signals. 

15 

BACKGROUND OF THE INVENTION 
Current pulse oximeters obtain two signals derived from the attenuation of red 
and infrared light signals as they are passed through a patient tissue site, typically a 
finger. A number of processing methods have been developed in the industry in both 

20 . time and frequency domains to obtain both pulse rate information and the oxygen content 
(Sp02) level of the arterial blood from the attenuated red and infrared light signals. The 
attenuated red and infrared signals show a pulsing waveform that is related to the heart 
rate of the patient. These time domain signals, usually after some bandpass filtering, are 
used for display of the pulse cycle and are known as plethysmographic signals. Prior 

25 techniques for pulse-rate estimation have mostly operated in the time domain and have 
used peak picking and analysis to derive a pulse rate. Time domain measures can 
respond quickly to pulse rate changes, but the presence of moderate motion and/or low 
amplitude pulses pose problems for accurate peak picking. Processing in the frequency 
or spectral domain has also been used and this requires a longer sample of the waveform 

30 to generate a pulse estimate. Also, identification of the predominant spectral peak 
produced by the pulse can be problematic in the presence of motion artifacts. 
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SUMMARY OF THE INVENTION 
Accordingly, the present invention provides for processing of plethysmographic 
signals via the cepstral domain to enhance the determination of patient physiological 
5 condition related information such as patient pulse rate and SP02 level information from 
plethysmographic signals, especially when motion artifacts are present in the 
plethysmographic signals. In accordance with the present invention, plethysmographic 
signals (e.g., attenuated red and infrared signals) are sampled and transformed into the 
cepstral domain, via, for example, a logarithmic like transform sandwiched between two 

10 forward Fourier transforms. Peaks in the cepstral domain are related primarily to the 
pulse rate. Identification of the pulse generated cepstral domain peak allows for pulse 
estimation in the presence of moderate motion artifacts. The cepstral information also 
allows for adaptive filtering of the input plethysmographic signals to remove noise and 
artifacts. The relative magnitudes of the cepstral peaks for both red and infrared signals 

15 in conjunction with an estimate of the DC levels of the red and infrared signals also 
allows for measurement of blood analyte (e.g., SP02) levels of the blood. 

According to one aspect of the present invention, a method of processing at least 
first and second time domain plethysmographic signals (e.g., red and infrared 
plethysmographic signals) obtained from a patient includes the steps of performing a 

20 Fourier transformation on the first time domain plethysmographic signal to transform the 
first plethysmographic signal into a first frequency domain plethysmographic signal and 
performing a Fourier transformation on the second time domain plethysmographic signal 
to transform the second plethysmographic signal into a second frequency domain 
plethysmographic signal. In this regard, the Fourier transformations may be fast Fourier 

25 transforms. A first power spectrum is computed from the first frequency domain 
plethysmographic signal and a second power spectrum is computed from the second 
frequency domain plethysmographic signal. A Fourier transformation is performed on 
the first power spectrum to transform the first power spectrum into a first cepstrum and a 
Fourier transformation is performed on the second power spectrum to transform the 

30 second power spectrum into a second cepstrum. In this regard, the Fourier 
transformations may be fast Fourier transforms. The first and second cepstrums are then 
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examined to obtain information therefrom relating to a physiological condition of the 
patient. 

The physiological condition of the patient may, for example, be the patient's pulse 
rate. In this regard, the first and second cepstrums may be examined to identify peaks in 
5 the first and second cepstrums associated with the pulse rate of the patient, and the pulse 
rate of the patient may be estimated based on the locations of the identified peaks in the 
first and second cepstrums. 

The physiological condition of the patient may also, for example, be the patient's 
SP02 level. In this regard, DC levels of the first and second power spectrums may be 

10 determined, AC levels of the first and second time domain plethysmographic signals may 
be determined from the identified peaks in the first and second cepstrums, and a value 
correlated with a blood analyte level (e.g., SP02 level) of the patient may be computed 
from the DC values of the first and second power spectrums and the AC levels of the first 
and second time domain plethysmographic signals. 

15 According to another aspect of the present invention, a method of determining a 

pulse rate of a patient from at least one time domain plethysmographic signal obtained 
from the patient includes the step of obtaining a time domain based estimate of the pulse 
rate of the patient from the time domain plethysmographic signal. The time domain 
plethysmographic signal is transformed to a spectral domain plethysmographic signal and 

20 a spectral domain based estimate of the pulse rate of the patient is obtained from the 
spectral domain plethysmographic signal. The spectral domain plethysmographic signal 
is transformed to a cepstral domain plethysmographic signal and a cepstral domain based 
estimate of the pulse rate of the patient is obtained from the cepstral domain 
plethysmographic signal. A best estimate of the pulse rate of the patient is then 

25 determined based on at least the time, spectral, and cepstral domain based estimates of 
the pulse rate of the patient. 

According to one more aspect of the present invention, a pulse oximeter includes 
first and second optical signal sources operable to emit optical signals characterized by 
first and second wavelengths (e.g., red and infrared), respectively. The pulse oximeter 

30 also includes a drive system, a detector, a digital sampler (e.g., an analog-to-digital 
converter), and a digital processor. The drive system is operable to cause operation of the 
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first and second optical signal sources such that each optical signal source emits first and 
second optical signals, respectively, in accordance with a multiplexing method. The 
detector is operable to receive the first and second optical signals after the first and 
second optical signals are attenuated by a patient tissue site of a patient. The detector is 
5 also operable to provide an analog detector output signal representative of the attenuated 
first and second optical signals. The digital sampler is operable to sample the analog 
detector output signal at a desired sampling rate and output a digital signal having a series 
of sample values representative of the attenuated first and second optical signals. The 
digital processor is enabled to demultiplex the series of sample values into first and 

10 second time domain plethysmographic signals, transform the first and second time 
domain plethysmographic signals into first and second spectral domain signals, transform 
the first and second spectral domain plethysmographic signals into first and second 
cepstral domain plethysmographic signals, and examine the first and second cepstral 
domain plethysmographic signals to obtain information therefrom relating to a 

15 physiological condition of the patient, such as the patient's pulse rate or SP02 level. 

According to a fiirther aspect of the present invention, a pulse arbitration method 
for use in determining a fiindamental pulse frequency (or pulse rate) of a patient firom 
multiple signal domains (e.g., time, energy, log, and cepstral) associated with at least one 
time domain plethysmographic signal obtained from the patient includes the step of 

20 transforming the time domain plethysmographic signal to a spectral domain 
plethysmographic signal. The spectral domain plethysmographic signal is transformed to 
a cepstral domain plethysmographic signal. The transformations to the spectral and 
cepstral domains may, for example, be accomplished via Fourier transformation 
operations. The spectral and cepstral domain plethysmographic signals are examined to 

25 identify corresponding spectral and cepstral domain plethysmographic signal peaks. The 
identified corresponding spectral and cepstral domain plethysmographic signal peaks are 
then used to select the fundamental pulse frequency from among a plurality of possible 
candidates for the fiindamental pulse frequency of the patient. Possible candidates for the 
fiindamental pulse frequency of the patient may, for example, be obtained from the time, 

30 spectral, and/or cepstral domain plethysmographic signals as well as from a filtered time 
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domain plethysmographic signal and/or a log scaled spectral domain plethysmographic 
signal. 

Cepstral domain processing of plethysmographic signals offers several advantages 
for pulse-rate identification. For example, the log-like transform acts to suppress weaker 
5 noise components making peak identification easier, and the cepstral peak is primarily 
generated via the harmonic components of the pulse so that noise energy surrounding this 
"fundamental" pulse frequency does not adversely effect pulse frequency identification. 

These and other aspects and advantages of the present invention will be apparent 
upon review of the following Detailed Description when taken in conjunction with the 
10 accompanying figures. 

DESCRIPTION OF THE DRAWINGS 
For a more complete understanding of the present invention and further 
advantages thereof, reference is now made to the following Detailed Description, taken in 
1 5 conjunction with the drawings, in which: 

FIG. 1 is a block diagram of one embodiment of a pulse oximeter in which a 
cepstral domain plethysmographic signal processing method in accordance with the 
present invention may be implemented; 

FIG. 2 is a block diagram showing one embodiment of a method for processing 
20 plethysmographic signals via the cepstral domain in accordance with the present 
invention; 

FIG. 3A is a plot showing typical red and infrared time domain plethysmographic 
input signals to be processed in accordance with the steps of FIG. 2; 

FIG. 3B is a plot showing the Spectrum and Cepstrum for the red 
25 plethysmographic input signal of FIG. 3 A after processing in accordance with the steps of 
FIG. 2; 

FIG. 3C is a plot showing the Spectrum and Cepstrum for the infrared 
plethysmographic input signal of FIG. 3 A after processing in accordance with the steps of 
FIG. 2; 
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FIG. 3D is a plot showing a typical infrared time domain plethysmographic input 
signal wherein the pulse oximeter probe is not transmitting properly through a patient 
tissue site (e.g., where the probe is removed from the patient's finger); 

FIG. 3E is a plot showing the Spectrum and Cepstrum for the infrared 
5 plethysmographic signal of FIG. 3D after processing in accordance with the steps of FIG. 
2; 

FIG.4 is a block diagram showing another embodiment of a method for 
processing plethysmographic signals via the cepstral domain in accordance with the 
present invention; 

10 FIG. 5 A is a plot showing typical red and infrared time domain plethysmographic 

input signals to be processed in accordance with the steps of FIG. 4; 

FIG. 5B is a plot showing differentiated waveforms obtained from the typical red 
and infrared time domain plethysmographic input signals shown in FIG. 5A; 

FIG. 5C is a plot showing red and infrared energy spectra corresponding to the 
15 typical red and infrared time domain plethysmographic input signals shown in FIG. 5 A; 
FIG. 5D is a plot showing red and infrared log spectra corresponding to the 
typical red and infrared time domain plethysmographic input signals shown in FIG. 5A; 

FIG. 5E is a plot showing red and infrared cepstrums corresponding to the typical 
red and infrared time domain plethysmographic input signals shown in FIG. 5A; 
20 FIG. 5F is a plot showing frequency domain filtered red and infrared 

plethysmographic waveforms corresponding to the typical red and infrared time domain 
plethysmographic input signals shown in FIG. 5A; 

FIG. 6A is a plot showing typical red and infrared time domain plethysmographic 
input signals to be processed in accordance with the steps of FIG. 4 that include motion 
25 induced noise components at a main motion frequency of about 200 bpm; 

FIG. 6B is a plot showing differentiated waveforms obtained from the typical red 
and infrared time domain plethysmographic input signals shown in FIG. 6A; 

FIG. 6C is a plot showing red and infrared energy spectra corresponding to the 
typical red and infrared time domain plethysmographic input signals shown in FIG. 6A; 
30 FIG. 6D is a plot showing red and infrared log spectra corresponding to the 

typical red and infrared time domain plethysmographic input signals shown in FIG. 6A; 
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FIG. 6E is a plot showing red and infrared cepstrums corresponding to the typical 
red and infrared time domain plethysmographic input signals shown in FIG. 6A; and 

FIG. 6F is a plot showing frequency domain filtered red and infrared 
plethysmographic waveforms corresponding to the typical red and infrared time domain 
5 plethysmographic input signals shown in FIG. 6A. 



DETAILED DESCRIPTION 
Referring now to Fig. 1, there is shown a block diagram of one embodiment of a 
pulse oximeter 10 in which a cepstral domain plethysmographic signal processing 

10 method in accordance with the present invention may be implemented. The pulse 

oximeter 10 is configured for use in determining the pulse rate of a patient as well as one 
or more blood analyte levels in the patient, such as an SP02 level. It should be 
appreciated that a cepstral domain plethysmographic signal processing method in 
accordance with the present invention may be implemented in pulse oximeters that are 

15 configured differently from the pulse oximeter depicted in FIG. 1 as well as in other 

environments wherein plethysmographic signals are processed in order to obtain desired 
information relating to patient physiological conditions from the plethysmographic 
signals. 

The pulse oximeter 10 includes a pair of optical signal sources 20a, 20b for 
20 emitting a corresponding pair of light signals 30a, 30b centered at different 

predetermined center wavelengths A i, /1 2 through a suitable tissue site of a patient and 
on to a detector 40 (e.g., a photo-sensitive diode). The optical signal sources 20a, 20b 
and detector 40 may be included in a positioning device 50, or probe, to facilitate 
alignment of the light signals 30a, 30b with the detector 40. For example, the positioning 
25 device 50 may be of clip-type or flexible strip configuration adapted for selective 

attachment to a suitable patient tissue site (e.g., a finger, an ear lobe, a foot, or the nose of 
the patient). The center wavelengths /1 1, A 2 required depend upon the blood analyte 
level to be determined. For example, in order to determine an SP02 level, A i may be in 
the Red wavelength range and X 2 may be in the infrared wavelength range. It should be 
30 appreciated that the pulse oximeter 10 may be readily implemented with more optical 
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signal sources (e.g., four) depending upon the number of different blood analyte levels to 
be measured. 

The optical signal sources 20a, 20b are activated by a corresponding plurality of 
drive signals 60a, 60b to emit the light signals 30a, 30b. The drive signals 60a, 60b are 
5 supplied to the optical signal sources 20a, 20b by a corresponding plurality of drive 
signal sources 70a, 70b. The drive signal sources 70a, 70b may be connected with a 
digital processor 80, which is driven with a clock signal 90 from a master clock 100. The 
digital processor 80 may be programmed to define modulation waveforms, or drive 
patterns, for each of the optical signal sources 20a, 20b. More particularly, the digital 

10 processor 80 may provide separate digital trigger signals 1 10a, 1 10b to the drive signal 
sources 70a-d, which in turn generate the drive signals 60a, 60b. In this regard, the 
digital trigger signals 1 10a, 1 10b may be configured to provide for multiplexing of the 
drive signals 60a, 60b, and in turn the light signals 30a, 30b, in accordance with a 
muhiplexing scheme (e.g., time division, frequency division, or code division 

15 multiplexing). 

The drive signal sources 70a, 70b, processor 80 and clock 100 may all be housed 
in a monitor unit 120. While the illustrated embodiment shows the optical signal sources 
20a, 20b physically interconnected with the positioning device 50 (e.g., mounted within 
the positioning device 50 or mounted within a connector end of a cable that is selectively 

20 connectable with the positioning device 50), it should be appreciated that the optical 
signal sources 20a, 20b may also be disposed within the monitor unit 120. In the latter 
case, the light signals 30a, 30b emitted from the optical signal sources 20a, 20b may be 
directed from the monitor unit 120 via one or more optical fibers to the positioning 
device 50 for transmission through the tissue site. Furthermore, the drive signal sources 

25 70a, 70b may comprise a single drive signal generator unit that supplies each of the drive 
signals 60a, 60b to the optical signal sources 20a, 20b. 

Transmitted light signals 130a, 130b (i.e., the portions of light signals 30a, 30b 
exiting the tissue) are detected by the detector 40. The detector 40 detects the intensities 
of the transmitted signals 130a, 130b and outputs a current signal 140 wherein the current 

30 level is indicative of the intensities of the transmitted signals 130a, 130b. As may be 
appreciated, the current signal 140 output by the detector 40 comprises a multiplexed 
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signal in the sense that it is a composite signal including information about the intensity 
of each of the transmitted signals 130a, 130b. Depending upon the nature of the drive 
signals 60a, 60b, the current signal 140 may, for example, be time division multiplexed, 
wavelength division multiplexed, or code division multiplexed. 
5 The current signal 140 is directed to an amplifier 150, which may be housed in the 

monitor unit 120 as is shown. As an alternative, the amplifier 150 may instead be 
included in a probe/cable unit that is selectively connectable with the monitor unit 120. 
The amplifier 150 converts the current signal 140 to a voltage signal 160 wherein a 
voltage level is indicative of the intensities of the transmitted signals 130a, 130b. The 

10 amplifier 150 may also be configured to filter the current signal 140 from the detector 40 
to reduce noise and aliasing. By way of example, the amplifier 150 may include a 
bandpass filter to attenuate signal components outside of a predetermined frequency 
range encompassing modulation frequencies of the drive signals 60a, 60b. 

Since the current signal 140 output by the detector 40 is a multiplexed signal, the 

15 voltage signal 160 is also a multiplexed signal, and thus, the voltage signal 160 must be 
demultiplexed in order to obtain signal- portions corresponding with the intensities of the 
transmitted light signals 130a, 130b. In this regard, the digital processor 80 may be 
provided with demodulation software for demultiplexing the voltage signal 160. In order 
for the digital processor 80 to demodulate the voltage signal 160, it must first be 

20 converted from analog to digital. Conversion of the analog voltage signal 160 is 

accomplished with an analog-to-digital (A/D) converter 170, which may also be included 
in the monitor unit 120. The A/D converter 170 receives the analog voltage signal 160 
from the amplifier 150, samples the voltage signal 160, and converts the samples into a 
series of digital words 180 (e.g., eight, sixteen or thirty-two bit words), wherein each 

25 digital word is representative of the level of the voltage signal 160 (and hence the 

intensities of the transmitted light signals 130a, 130b) at a particular sample instance. In 
this regard, the A/D converter 170 should provide for sampling of the voltage signal 160 
at a rate sufficient to provide for accurate tracking of the shape of the various signal 
portions comprising the analog voltage signal 160 being converted. For example, the 

30 A/D converter 170 may provide for a sampling frequency at least twice the frequency of 
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the highest frequency drive signal 60a, 60b, and typically at an even greater sampling rate 
in order to more accurately represent the analog voltage signal. 

The series of digital words 180 is provided by the A/D converter 170 to the 
processor 80 to be demultiplexed. More particularly, the processor 80 may periodically 
5 send an interrupt signal 190 (e.g., once per every eight, sixteen or thirty- two clock cycles) 
to the A/D converter 170 that causes the A/D converter 170 to transmit one digital word 
180 to the processor 80. The demodulation software may then demultiplex the series of 
digital words 180 in accordance with an appropriate method (e.g., time, wavelength, or 
code) to obtain digital signal portions indicative of the intensities of each of the 

10 transmitted light signals 130a, 130b. In this regard, the demuhiplexed digital signal 
portions comprise time domain plethysmographic signals corresponding to the center 
wavelengths /I i, /I2 (e.g., red and infrared) of the optical signal sources 20a, 20b. The 
red and infrared time domain plethysmographic signals may then be processed by the 
processor 80 to obtain desired patient physiological condition related information 

15 therefrom such as the patient's pulse rate and SP02 level. 

Referring now to FIG. 2 there is shown a block diagram illustrating one 
embodiment of a method (200) for processing the red and infrared time domain 
plethysmographic signals via the cepstral domain to obtain desired information relating to 
patient physiological conditions such as patient pulse rate and blood analyte level (e.g., 

20 SP02) information. The cepstral domain plethysmographic signal processing method 
(200) begins with obtaining (210) two digitized time domain plethysmographic signals 
such as red and infrared plethysmographic signals. In this regard, typical red and infrared 
time domain plethysmographic signals that have been sampled at 50 Hz are shown in 
FIG. 3A. The cepstral domain processing method (200) is particularly suited for 

25 implementation in software executable by the digital processor 80 of a pulse oximeter 10 
such as described above in connection with FIG. 1. In other embodiments, the cepstral 
domain processing method (200) may be configured for processing non-digitized 
plethysmographic signals and may be implemented in appropriate hardware components. 
Furthermore, the cepstral domain processing method (200) may be configured for 

30 simultaneously processing more than two plethysmographic signals. 
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A suitable smoothing window function (e.g., Hanning, Hamming, Kaiser) is 
applied (220) to the digitized time domain plethysmographic signals to smooth the 
signals. Smoothing the digitized time domain plethysmographic signals achieves 
improved frequency estimation. After the signals are smoothed, a first Fourier 
5 transformation operation is performed (230) on the signals to transform the red and 
infrared plethysmographic signals from the time domain to the frequency domain. Since 
there are two primary signals (the red and infrared inputs), it is convenient to perform the 
first Fourier transformation of the signals in parallel using a complex Fast Fourier 
Transform (FFT) procedure. If desired, the results of the FFT calculations may be 

10 appropriately scaled (e.g., by dividing by the number of points used in the FFT 
calculations) to help prevent floating point overflow errors in subsequent computations. 
After the first stage FFT is performed, respective power spectrums are computed (240) 
from the frequency domain red and infrared plethysmographic signals. In this regard, the 
power spectrums may be computed (240) by squaring and summing the appropriate real 

15 and imaginary frequency components of the red and infrared frequency domain 
plethysmographic signals. Power spectrums of the typical red and infrared 
plethysmographic signals after the first stage FFT are shown in FIGS. 3B and 3C, 
respectively. 

After the power spectrums are computed, a log-like or companding function is 
20 applied (250) to the red and infrared power spectrums. Application of the log-like or 
companding function suppresses smaller noise components and emphasizes the 
prominent harmonics so that periodicity in the spectrum is more easily extracted. A 
second Fourier transformation operation is then performed (260) on the log transformed 
power spectrums to transform the signals to the cepstral domain. In this regard, it is 
25 convenient to perform the second-stage Fourier transformation of the log scaled power 
spectrums in parallel using a complex Fast Fourier Transform (FFT) procedure. If 
desired, the results of the second-stage FFT calculations may be appropriately scaled in a 
manner similar to scaling done on the results of the first-stage FFT calculations. The 
cepstrums of the typical red and infrared plethysmographic signals obtained after the 
30 second stage FFT are also shown in FIGS. 3B and 3C, respectively. 



11 



Attorney Docket No. 41942-05553 

Once the red and infrared cepstrums are obtained, the separate red and infrared 
cepstrums are then examined (270) for peaks associated with the pulse rate of the patient. 
In this regard, the most prominent (i.e., largest amplitude) peak in each cepstrum may be 
identified. The location of the most prominent peak in each cepstrum provides an 
5 indication of the fundamental frequency of the plethysmographic waveform from which 
the cepstrum is obtained. Since the fundamental frequency of a plethysmographic 
waveform is proportional to the patient's pulse rate, the pulse rate of the patient may be 
estimated (280) from one or both of the cepstrums. For example, the most prominent 
peak in the red cepstrum of FIG. 3B occurs at around the 20th bin of the FFT spectrum 

10 corresponding to a cepstral based pulse rate estimate of approximately 65 beats-per- 
minute. It should be noted that this estimate differs slightly from a conventional time 
domain based estimate obtained from the time domain red plethysmographic waveform 
shown in FIG. 3 A of 61 beats-per-minute. Pulse-rate estimates may be obtained from 
both the red and infrared cepstrums and the separate estimates may be correlated with one 

15 another in order to obtain a single estimate of the patient's pulse rate. Further, while it is 
possible to estimate the patient's pulse rate based only on information from one or both 
of the cepstrums, a time domain based estimate of the patient's pulse rate may also be 
used for initial identification purposes and to support subsequent tracking of the cepstral 
peak (Quefrency) associated with the pulse rate. 

20 In some cases, there may not be a prominent peak in one or both of the cepstrums. 

For example, FIG. 3D shows an infrared time domain plethysmographic signal typical of 
the situation where there is no physiological signal condition (e.g., where the 
plethysmographic probe has been removed from the patient's finger), and FIG. 3E shows 
the infrared power spectrum and cepstrum obtained for the infrared time domain 

25 plethysmographic signal of FIG. 3D. While the power spectrum of FIG. 3E differs 
somewhat from a power spectrum that is typical of a patient physiological signal 
condition such as the power spectrums shown in FIGS. 3B and 3C, the lack of a patient 
physiological signal condition is particularly apparent from examination of the cepstrum 
since there is no prominent peak present in the cepstrum of FIG. 3E as compared with the 

30 quite prominent cepstral peaks in FIGS. 3B and 3C. 
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In addition to examining the cepstrums for peaks associated with patient pulse 
rate, in step (270) the red and infrared cepstrums may be examined for peaks associated 
with motion artifacts. Typically, peaks in the red and infrared cepstrums that are 
associated motion artifacts will be less prominent than the peaks associated with the 
5 patient pulse rate. The location(s) of less prominent peaks in each cepstrum provide an 
indication regarding motion artifacts present in the plethysmographic waveform from 
which the cepstrum is obtained, and based on this information the frequencies of motion 
artifacts present in the red and infrared plethysmographic signals may be estimated (290). 
Once an estimate of the pulse rate is obtained, the pulse rate information may be 

10 used to construct a filter to remove noise and motion artifacts from the input red and 
infrared signals. This may be done via an adaptive bandpass filter applied in the time 
domain to the red and infrared signals where the cut off frequencies are determined by 
the pulse frequency which is identified in the cepstral domain. Alternatively, as is shown 
in the embodiment of FIG. 2, the frequency domain red and infrared plethysmographic 

15 signals may be filtered (300) in the frequency domain after the first stage FFT with a 
frequency domain fiher constructed using the pulse frequency information obtained from 
the cepstral domain. An inverse fast Fourier transform (IFFT) operation may be 
performed (310) on the filtered frequency domain signals to obtain filtered time domain 
red and infrared plethysmographic signals for use in subsequent measures such as a 

20 regression based SP02 estimation which uses the time domain version of the red and 
infrared inputs signals. Noise removal from the red and infrared signals improves 
subsequent measures such as regression based SP02 estimation. 

Additionally, the information in both the spectral and cepstral domains may be 
used to derive an SP02 measure. The overall DC levels of the red and infrared 

25 plethysmographic signals can be determined from the first stage spectrums and the 
relative magnitudes of the cepstral peaks corresponding to the pulse rate frequency may 
be used to obtain a measure of the AC levels of the red and infrared plethysmographic 
signals. In this regard, the following computation may be utilized: 

30 ^ ' = AC(cepstral'red) / DCfspectral-red) / ACfcepstral-IR) / DCfspectral-IR) 
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or, expressed in an alternative manner: 

R' ^ AC(cepstral-red) / DC(spectral'red) * DC(spectral-IR) / AC(cepstral-IR) 

5 where ACfcepstral-red) is the AC level of the red plethysmographic signal obtained from 
the red cepstrum, DC(spectral-red) is the DC level of the red plethysmographic signal 
obtained from the red spectrum, AC(cepstral-IR) is the AC level of the infrared 
plethysmographic signal obtained from the infrared cepstrum, and DC(spectral-IR) is the 
DC level of the infrared plethysmographic signal obtained from the infrared spectrum. 

10 The derived measure R' may then be used to estimate (320) the patient's SP02 level in a 
manner similar to known regression techniques where AC and DC estimates are obtained 
from the time domain red and infrared signals. An example of such a known regression 
technique is described in United States Patent No. 5,934,277 entitled "SYSTEM FOR 
PULSE OXIMETRY SP02 DETERMINATION", the entire disclosure of which is 

1 5 incorporated herein. 

Referring now to FIG. 4 there is shown a block diagram illustrating another 
embodiment of a method (400) for processing the red and infrared time domain 
plethysmographic signals via the cepstral domain to obtain desired information relating to 
patient physiological conditions such as patient pulse rate and blood analyte level (e.g., 

20 SP02) information. The cepstral domain plethysmographic signal processing method 
(400) shown in FIG. 4 proceeds in a manner similar to the method (200) shown in FIG. 2. 
In this regard, two continuous time domain plethysmographic signals such as red and 
infrared plethysmographic signals are digitized (410) by sampling the signals at a suitable 
frequency. Typical red and infrared time domain plethysmographic signals that have 

25 been sampled at 50 Hz are shown in FIGS. 5 A and 6A, with the signals of FIG. 6A 
including motion artifacts. As with the method (200) of FIG. 2, the cepstral domain 
processing method (400) is particularly suited for implementation in software executable 
by the digital processor 80 of a pulse oximeter 10 such as described above in connection 
with FIG. 1, and in other embodiments, the cepstral domain processing method (400) 

30 may be configured for processing non-digitized plethysmographic signals and may be 
implemented in appropriate hardware components. Furthermore, the cepstral domain 
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I processing method (400) may be configured for simultaneously processing more than two 

' plethysmographic signals. 

The digitized time domain red and infrared plethysmographic signals are 
smoothed (420) via a suitable smoothing window (e.g. Banning, Hamming, or Kaiser) 
5 and are then processed in parallel via a complex FFT (430). The output from the first 
stage FFT is then decoded and the separate red and infrared energy spectra and log power 
spectra are computed and stored (440, 450). Plots of red and infrared energy spectra and 
log spectra obtained for the red and infrared signals of FIG. 5 A and 6A are shown in 
FIGS. 5C and 5D, respectively, and in FIGS. 6C and 6D, respectively. A second stage 
10 FFT (460) is then applied to the log power spectra to obtain red and infrared cepstra (470) 
therefrom. If desired, the results of the first and second stage FFT calculations may be 
scaled to help prevent floating point errors in subsequent computations. Plots of the red 
and infrared cepstra obtained for the red and infrared signals of FIG. 5 A and 6A are 
shown in FIG. 5E and 6E. Peaks in the cepstra (which has the dimension of Quefrency) 
15 are examined (480) and transformed to provide an estimate of pulse frequency. 

The cepstral based pulse rate estimate is provided to a pulse arbitration module 
(490). The pulse arbitration module (490) also receives estimates of the patient's pulse 
rate based on examination of peaks in the energy spectra and log power spectra. 
Additionally, a time-domain pulse rate estimate is extracted (500) from the digitized time 
20 domain red and infrared plethysmographic signals via a conventional technique such as 
differentiation , thresholding and picking the most commonly found interval. Plots of 
the differentiated waveforms obtained from the time domain red and infrared 
plethysmographic signals of FIGS. 5 A and 6A are shown in FIGS. 5B and 6B. The time 
domain based pulse rate estimate is also provided as an input to the pulse arbitration 
25 module (490). 

Information relating to the peaks of the energy spectra and the cepstra, are input to 
a motion classification and motion strength estimation module (510). The motion 
classification and motion strength estimation module (510) uses both the amplitude, 
relative position and spacing of the respective peaks in the red and infrared energy 
30 spectra and cepstra to make motion classification and strength judgments. A simple 

measure classification and motion estimation can be derived by the number and spacing 
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of cepstral peaks. In this regard, a relatively clean plethysmographic signal will typically 
produce one major cepstral peak. As the number and size of the cepstral peaks increases, 
sizable motion components can be inferred. Information from the motion classification 
and motion strength estimation module (510) is input to both an adaptive filter module 
5 (520) and the pulse arbiter module (490). 

The adaptive filter module (520) uses estimates of the pulse frequency and the 
frequency distribution of the motion noise components (if present) to control filtering in 
the frequency domain in order to improve the signal to noise ratio of the pulse 
fundamental frequency components and/or its harmonics. In this regard, the red and 

10 infrared frequency domain plethysmographic signals obtained after the first stage FFT 
(430) signals are filtered (530) to produce filtered frequency domain red and infrared 
plethysmographic signals. Plots of the filtered frequency domain red and infrared 
plethysmographic signals corresponding to the time domain red and infrared 
plethysmographic signals of FIGS. 5A and 6A are shown in FIGS. 5F and 6F. A 

15 number of different types of filters may be implemented including both finite impulse 
response (FIR) and infinite impulse response (IIR) filters. One disadvantage of spectral 
methods is that they are not suited for tracking rapid changes in the input signal. 
However in the present method (400) the spectral information is used to control an 
adaptive filter. By using time domain pulse measurement techniques on the output signal 

20 from this filter, the ability to track reasonably fast changes is achieved. 

An inverse FFT operation (540) is performed to obtain filtered time domain red 
and infrared plethysmographic signals, and an overlap and add operation (550) is 
performed to reconstruct the plethysmographic signals minus the DC components and 
with reduced motion components. Following the overlap and add operation (550), the 

25 energy content for both the red and infrared filtered signals is then obtained (560) via, for 
example, a root-mean-square (rms) measure. This provides an estimate of the AC red 
and infrared levels. Although not shown in FIG. 4, it is also possible to obtain an 
estimate of the red and infrared AC levels via the cepstral domain. The main peak 
location of the red and infrared cepstra can be translated to a frequency value and the 

30 value of the energy for that frequency and its harmonics can be obtained (i.e., integrated) 
by referring to the stored energy spectrum for the red and infrared signals. It is also 
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feasible to use the relative amplitudes of the red and infrared cepstral peaks to derive an 
AC estimate. Following the overlap and add operation (550), another conventional time 
domain based pulse estimation is also performed (570) on the filtered red and infrared 
signals and this estimate is also sent to the pulse arbiter module (490). 
5 The pulse arbiter (490) uses the various time domain, filtered time domain, energy 

spectra, log power spectra and cepstral based pulse estimates and the motion strength and 
classification to provide an overall best estimate (580) of the patient's pulse rate. In this 
regard, for a range of motions the location of the major cepstral peak suffices as a good 
estimate of pulse frequency. However for large motion amplitudes and motion that 

10 produces waveforms similar to those of red and infrared plethysmographic signals it is 
necessary to examine a number of parameters to resolve competing estimates. More 
particularly, the pulse arbitration module (490) examines the correlation between the time 
domain (both filtered and unfiltered), spectral domain (both energy and log power) and 
cepstral domain based pulse estimates and uses the motion estimation derived from the 

15 cepstrum in the motion classification and motion strength estimation module (510) to 
weight the respective pulse rate estimates. If significant motion is present then cepstral 
information can be used to resolve between competing spectral pulse candidates. In this 
regard, the typical pulse waveform which is 'sawtoothed shaped' would result in a main 
fundamental spectral peak with usually at least two visible harmonic peaks. The resulting 

20 cepstra would be one main peak associated with the fundamental frequency. Therefore 
spectral candidates with no corresponding cepstral peak can be eliminated. Further, in 
cases where there is competing noise around the fundamental frequency peak, a cepstral 
peak can be confirmed by examining the energy or log spectra for a fundamental 
frequency peak and related harmonics. In addition to the previously described pulse 

25 arbitration process (490), it would also be feasible to employ a neural-net for the pulse 
arbitration process (490). 

Another strategy that may be employed in the pulse arbiter module (490) is to 
relate the cepstral peak to a region or channel in the energy spectrum and to obtain an AC 
value and then derive a SP02 estimate. This SP02 estimate can be referred to another 

30 SP02 estimate derived from the mean energy over the allowable pulse range (e.g., 30- 
350 bpm). A valid cepstral candidate will generate a similar track of SP02 over time as 
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the estimate derived from mean energy. This information can also be used to resolve 
amongst competing cepstral candidates for the one related to the pulse frequency. 

In addition to obtaining an overall best estimate (580) of the patient's pulse rate, 
the plethysmographic signal processing method (400) of FIG. 4 also derives an estimate 
5 of the patient's SP02 level. The energy content of the time domain red and infrared 
plethysmographic signals is obtained (590) via, for example a root mean square (rms) 
transform. This provides an estimate of the red and infrared DC levels. The red and 
infrared DC levels (590) and AC levels (560) are provided to an SP02 module (600). As 
discussed in more detail above in connection with step (310) of the method (200) of FIG. 

10 2, the SP02 module (600) uses the red and infrared DC and AC levels to derive a 
measure that can be correlated with the patient's SP02 level in a manner similar to 
conventional regression based techniques. 

The cepstral domain plethysmographic signal processing method (400) of FIG. 4 
also provides for obtaining an enhanced perfusion index (PI) measure when motion 

15 artifacts are present in the red and infrared time domain plethysmographic signals as 

compared to known time domain based perfusion index measures. The perfusion index 
is a measure of relative perfusion in the patient tissue site and is indicative of pulse 
strength. A time-domain based perfusion index measure may be obtained by, for 
example, calculating normalized plethysmographic signal amplitudes for the red and 

20 infrared time domain plethysmographic signals by summing the normalized delta 

amplitudes covering the rising portion of one cycle of the pulse waveform. This value 
can be termed Snda. In this regard, the perfusion index may be calculated from the red 
and infrared Snda values in accordance with the following expression: 



25 PI = (Snda (red) * 0.0563 + Snda (infrared) * 0.3 103) * Scaling Factor 

Further detail regarding such a known time domain based method for obtaining a 
perfusion index measure is described in United States Patent No. 5,766,127 entitled 
"METHOD AND APPARATUS FOR IMPROVED PHOTOPLETHYSMOGRAPHIC 
30 PERFUSION-INDEX MONITORING", the entire disclosure of which is incorporated 
herein. 
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However it is also possible to obtain a measure of the red and infrared 
plethysmographic signal amplitudes from their respective energy spectrums when the 
frequency components present in the energy spectrums due to the pulse signal can be 
identified via processing of the red and infrared cepstrums. In the regard, the 
5 plethysmographic signal processing method (400) may incorporate a perfusion index 
estimator step (610) wherein the red and infrared cepstrums obtained in step (470) are 
used to identify the frequency components present in the red and infrared energy 
spectrums obtained in step (440) that are associated with the pulse rate of the patient (i.e. 
the fundamental pulse frequency and its harmonics). The perfusion index estimator 
10 module (610) computes normalized amplitudes for the identified red and infrared spectral 
peaks. A perfusion index value (620) may then be computed from the normalized 
amplitudes of the identified red and infrared spectral peaks in accordance with, for 
example, the following expression: 



15 PI = (ESamp (red) * 0.0563 + ESamp(infrared) + 0.3 103) * ESscaling 

where ESamp(red) and ESamp(infrared) are the normalized amplitudes derived from the 
identified spectral peaks in the red and infrared energy spectrums and ESscaling is a 
scaling factor adjusted to give the spectral PI measure an equivalent value to the time 

20 domain PI measure. Because the spectral PI measure uses normalized amplitudes of the 
identified peaks in the red and infrared spectrums associated with the fundamental pulse 
frequency, the spectral PI measure is less susceptible to corruption by motion artifacts 
present in the time domain plethysmographic signals since peaks associated with motion 
artifacts will be ignored when identifying the fundamental pulse frequency peaks using 

25 the cepstrums. 

Where desired, the spectral based PI measure may be correlated with the time 
domain based PI measure to provide a single PI measure. Further the spectral based PI 
measure provides information that can be used in tracking and identification of the 
fundamental pulse frequency by the pulse arbitration module (490). In this respect a 
30 spectral PI measure may be calculated for each spectral candidate and these estimates can 
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be used in a scoring and arbitration scheme to track and resolve the correct (pulse 
produced) fundamental pulse frequency candidate. 

While various embodiments of the present invention have been described in 
detail, further modifications and adaptations of the invention may occur to those skilled 
5 in the art. However, it is to be expressly understood that such modifications and 
adaptations are within the spirit and scope of the present invention. 
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